UNCLASSIFIED 


_ AD  NUMBER _ 

ADB383242 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  DoD  and  DoD 
contractors  only;  Administrative/Operational 
Use;  JUN  2012.  Other  requests  shall  be  referred 
to  Naval  Air  Systems  Command,  Code  4.3.1,  48110 
Shaw  Road,  Patuxent  River,  MD  20670-1906. 


_ AUTHORITY 

NAWCADPAX  errata  dtd  2  Jul  2013 


THIS  PAGE  IS  UNCLASSIFIED 


NAVAL  AIR  WARFARE  CENTER  AIRCRAFT  DIVISION 
PATUXENT  RIVER,  MARYLAND 


TECHNICAL  REPORT 


REPORT  NO:  NAWCADPAX/TR-2012/195 


POTENTIAL  CRASH  LOCATION  (PCL)  MODEL 

by 

Natasha  Bradley 
Dr.  David  Burke 

6  June  2012 


Distribution  authorized  to  DOD  and  U.S.  DOD  contractors  only;  Administrative 
or  Operational  Use;  June  2012.  Other  requests  shall  be  referred  to  the  Naval  Air 
Systems  Command  (Code  4.3.1),  48110  Shaw  Road,  Patuxent  River,  Maryland 
20670-1906. 

DESTRUCTION  NOTICE  -  Destroy  by  any  method  that  will  prevent  disclosure 
of  contents  or  reconstruction  of  the  document. 


UNCLASSIFIED 


DEPARTMENT  OF  THE  NAVY 
NAVAL  AIR  WARFARE  CENTER  AIRCRAFT  DIVISION 
PATUXENT  RIVER,  MARYLAND 


NAWCADPAX/TR-20 12/195 
6  June  2012 


POTENTIAL  CRASH  LOCATION  (PCL)  MODEL 
by 

Natasha  Bradley 
Dr.  David  Burke 


RELEASED  BY: 


ROLAND  COCHRAN  /  AIR-4.3. 1  /  DATE 

Systems  Engineering  Division 

Naval  Air  Warfare  Center  Aircraft  Division 


NAWCADPAX/TR-20 12/195 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering 
and  maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188), 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any 
penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE 
ADDRESS. 


1 .  REPORT  DATE  2.  REPORT  TYPE  3.  DATES  COVERED 

6  June  20 1 2  Technical  Report  June  -  November  2011 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 


Potential  Crash  Location  (PCL)  Model 


5b.  GRANT  NUMBER 


6.  AUTHOR!  S) 


5c.  PROGRAM  ELEMENT  NUMBER 
PE064400D8Z 


5d.  PROJECT  NUMBER 


Natasha  Bradley 
Dr.  David  Burke 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 


8.  PERFORMING  ORGANIZATION  REPORT  NUMBER 


Naval  Air  Warfare  Center  Aircraft  Division  (AIR-4.3.1) 

48110  Shaw  Road 

Patuxent  River.  Maryland  20670-1 161 _ 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

OUSD(AT&L)/S&TS-Unmanned  Warfare 
Room  3B938 
3090  Defense  Pentagon 
Washington,  DC  20301-3090 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 


NAWCADPAX/TR-2012/195 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT  NUMBER! S) 


Distribution  authorized  to  DOD  and  U.S.  DOD  contractors  only;  Administrative  or  Operational  Use;  June  2012.  Other  requests  shall  be  referred 
to  the  Naval  Air  Systems  Command  (Code  4.3.1),  48110  Shaw  Road,  Patuxent  River,  Maryland  20670-1906 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

As  the  use  of  Unmanned  Aerial  Systems  (UAS)  grows,  there  is  a  need  to  integrate  UAS  into  the  National  Airspace  System.  One  challenge  that 
has  been  encountered  is  developing  a  comprehensive  risk  analysis  approach  to  model  the  risk  posed  by  UAS  operations  to  the  general  public. 
One  component  of  this  Target  Level  of  Safety  approach  is  a  model  for  determining  crash  location.  There  are  multiple  approaches  for  modeling 
the  potential  crash  locations  of  an  UAS  after  a  catastrophic  failure.  This  report  examines  two  potential  models  for  fixed  wing  based  on  previous 
research  into  UAS  risk  analysis:  the  Clothier  model  (reference  1 )  and  the  Sensis  model  (reference  2)  as  well  as  two  models  for  rotary  wing 
aircraft. 


15.  SUBJECT  TERMS 


Potential  Crash  Location  (PCL);  Unmanned  Aerial  System  (UAS);  Clothier  model;  Sensis  model 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 


Unclassified  Unclassified 


Unclassified 


17.  LIMITATION  18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

OF  ABSTRACT  OF  PAGES  Dr.  David  Burke _ 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

SAR  22  (301)342-2185 


1 


NAWCADPAX/TR-20 12/195 


SUMMARY 

As  the  use  of  Unmanned  Aerial  Systems  (UAS)  grows,  there  is  a  need  to  integrate  UAS  into  the 
National  Airspace  System.  One  of  the  critical  requirements  for  expanding  the  operational  area  of 
UAS  is  to  understand  the  risk  to  uninvolved  third  parties  on  the  ground  posed  by  the  crash  of  a 
UAS.  In  order  to  address  this  issue  Office  of  Secretary  of  Defense  (OSD),  Strategic  and  Tactical 
Systems  -  Unmanned  Warfare  Office,  Air  Worthiness  IPT  has  sponsored  the  Target  Level  of 
Safety  to  Third  Parties  program.  The  objective  of  this  program  is  to  define  a  consistent 
calculation  method  to  detennine  the  relationship  between  UAS  reliability,  potential  to  cause 
damage  and  where  it  flies.  Naval  Air  Systems  Command,  AIR  4.3.1  has  lead  the  effort  to 
develop  this  methodology.  The  Target  Level  of  Safety  (TLS)  Program  includes  5  modules; 
Casualty  Expectation,  Probability  of  Loss  of  Aircraft,  Potential  Crash  Location,  Lethal  Crash 
Area  and  Population  Density.  This  paper  focuses  on  methods  to  define  the  Potential  Crash 
Location  of  a  UAS  given  a  catastrophic  failure. 


There  are  multiple  approaches  for  modeling  the  potential  crash  locations  of  an  UAS  after  a 
catastrophic  failure.  This  report  examines  two  potential  models  for  fixed  wing  aircraft  based  on 
previous  research  into  UAS  risk  analysis:  the  Clothier  model  (reference  1)  and  the  Sensis  model 
(reference  2).  Clothier  uses  a  high  fidelity  dynamic  model  of  a  generic  aircraft,  while  Sensis  Inc. 
demonstrates  a  worst-case/best-case  approach  to  defining  the  outer  limits.  This  report  also 
discusses  two  different  approaches  to  modeling  rotor  craft  UAS  crash  locations. 
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INTRODUCTION 


BACKGROUND 


As  the  use  of  Unmanned  Aerial  Systems  (UAS)  grows,  there  is  a  need  to  integrate  UAS  into  the 
National  Airspace  System.  One  challenge  that  has  been  encountered  is  developing  a 
comprehensive  risk  analysis  approach  to  model  the  risk  posed  by  UAS  operations  to  the  general 
public.  The  current  effort  is  examining  a  Target  Level  of  Safety  (TLS)  approach  to  modeling  this 
risk  to  third  parties,  those  not  involved  with  UAS  operations.  One  component  of  this  TLS 
approach  is  a  model  for  detennining  crash  location.  The  Potential  Crash  Location  (PCL)  model 
uses  some  simplified  equations  to  model  the  situation  where  an  UAS  has  a  catastrophic  failure  at 
some  altitude.  It  provides  the  footprint  within  which  the  Unmanned  Aircraft  (UA)  is  most  likely 
to  crash.  The  model  takes  into  account  both  fixed  wing  and  rotor  wing  aircraft  of  various 
configurations  and  sizes.  The  Lethal  Crash  Area  (LCA)  model  provides  a  method  for  calculating 
how  large  the  actual  lethal  area  at  the  site  of  impact  would  be.  The  LCA  model  is  described  in  a 
separate  report  (reference  3). 

There  are  several  published  approaches  for  modeling  the  potential  crash  locations  of  an  UAS 
after  a  catastrophic  failure.  This  report  examines  two  potential  models  for  fixed  wing  aircraft 
based  on  previous  research  into  UAS  risk  analysis:  the  Clothier  model  (reference  1)  and  a  Sensis 
model  (reference  2).  Clothier  uses  a  high  fidelity  dynamic  model  of  a  generic  aircraft,  while 
Sensis  Inc.  demonstrates  a  worst-case/best-case  approach  to  defining  the  outer  limits.  Both  of 
these  models  are  confined  to  fixed  wing  aircraft.  This  report  also  discusses  two  different 
approaches  to  modeling  rotor  craft  UAS  crash  locations. 

CLOTHIER  MODEL 


Clothier  (reference  1)  uses  a  6-degree  of  freedom  (6-DOF)  model  to  find  the  footprint  boundaries 
for  gliding  decent  using  differing  aircraft  parameters.  Clothier  uses  a  simple  6-DOF  model  to 
find  the  maximum  extent  of  the  crash  footprint  boundaries.  Clothier  first  makes  some 
assumptions  to  simplify  the  equations  of  motion.  Since  the  aircraft  is  a  UAS,  it  is  safe  to  assume 
that  it  has  coordinated  flight.  Clothier  also  assumes  that  the  transitions  between  trim  states  are 
instantaneous.  Therefore,  the  trajectories  this  model  generates  represent  an  upper  bound  on  the 
actual  trajectories  that  can  be  achieved.  The  model  also  assumes  simple  wing  geometry,  constant 
descent  velocity,  balanced  flight  (lift  force  equals  weight),  and  no  wind  effects.  The  outer  edge 
of  the  footprint  is  obtained  by  simulating  a  glide  descent  over  all  heading  angles  where  a  heading 
of  0  deg  corresponds  to  an  impact  directly  in  front  of  the  aircraft.  Note  however  that  a  heading  of 
180  deg  does  not  correspond  to  an  impact  point  directly  behind  the  aircraft  due  to  a  non- 
negligible  turn  radius.  Therefore,  Clothier  assumes  that  the  descent  trajectory  comprises  a  turn 
followed  by  a  straight  line  glide.  In  regards  to  this  turn,  Clothier  also  assumed  that  there  is  a 
constant  flight  path  angle  (FPA)  and  the  bank  angle  is  45  deg.  He  verified  the  assumption  of  a  45 
deg  bank  angle  by  checking  the  displacement  of  the  ground  impact  point  given  identical  initial 
conditions;  a  45  deg  bank  angle  gives  an  approximately  equal  displacement.  Clothier  also  made 
some  assumptions  that  made  the  mathematics  manageable  such  as  rigid  body  dynamics,  small 
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angle  of  attack  (AOA),  small  thrust  angle,  and  small  FPA.  A  key  assumption  is  that  the  entire 
flight  will  occur  at  sea  level  standard  day.  This  assumption  greatly  simplifies  the  equations  since 
varying  atmospheric  density  would  require  flight  by  flight  model  rather  than  a  generic  model. 
This  assumption  is  acceptable  for  this  tool  since  most  of  the  UAS  will  be  flying  at  fairly  low 
altitudes  and  any  UAS  up  at  higher  altitudes  would  have  a  shorter  glide  distance  due  to  lower  air 
density. 


In  order  to  find  the  maximum  extent  of  the  footprint  boundary,  the  FPA  is  needed.  That  can  be 
found  using  the  above  assumptions  with  the  differential  equations  for  velocity  and  FPA  (y).  The 
differential  equations  are  (reference  1): 


Y  = 


\t\ 

_ g_ 

Wh  x 


[||fi,||  x  cos(cr  +  e)  —  ||DW||  +  \\Wh\\  x  siny]  =  0 
Tb ||  x  sin(a  +  s)  +  ||ZW||)  x  cos fi  —  \\Wh\\  x  cosy]  =  0 


When  simplified  using  the  assumption  of  a  gliding  decent,  the  differential  velocity  would 
become: 


=  0 


By  substituting  the  drag  force  equation,  the  FPA  would  be: 

1  „  „  1 1  1 1 2 

y  = 


2pSCD  X  Vw 


m 


h\ 


The  AOA  given  y,  Vw  and  simplifying  differential  equation  for  y,  shown  above,  would  become: 

\\LW\\  cos  [a  -  \\Wh\\  =  0 


The  Coefficient  of  Lift  (Cl)  and  Coefficient  of  Drag  (Cd)  can  be  simplified  to  the  equations 
below: 

—  CiUq  +  Ciacc 

Cd  =  Cdo+7Xre 


The  6-DOF  model  parameters  can  be  evaluated  using  the  above  equations. 


It  is  important  to  note  that  the  minimum  value  for  y  is  achieved  when  the  aircraft  flies  at  its 
maximum  lift  to  drag  ratio.  That  ratio  is  achieved  at  the  minimum  drag  velocity: 
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^Min  Drag  lr  II 


2  W 


1  0.25 

t  ECnJ 


pS  \tcAR 

In  flying  this  maximum  lift  to  drag  ratio  profile,  the  maximum  gliding  range  is: 


cW  = 


H 


tan (ymin) 


H  C, 

- ~  H7^ 

Ymin  ld 


MAX 


Next,  Clothier  evaluates  his  maximum  boundary  results  by  plugging  in  values  from  different 
fixed  wing  UAV’s  into  the  6-DOF  model.  Clothier  found  that  footprint  area  grows  uniformly  as 
the  height  above  ground  level  increases.  At  low  heights,  the  footprint  is  more  like  a  pie  slice  and 
as  the  height  increases  the  footprint  looks  more  like  an  ellipse.  The  footprint  shape  and  area  also 
varies  with  different  aircraft,  mainly  due  to  the  difference  in  glide  ratios.  The  better  the  aircraft’s 
glide  ratio,  the  larger  the  footprint  area.  Initial  velocity  is  another  factor  that  determines  shape 
since  the  higher  the  velocity,  the  less  circular  the  footprint  becomes  due  to  the  limitation  on  bank 
angle.  Glide  velocity  has  the  opposite  effect  of  height  such  that  the  higher  the  glide  velocity  the 
more  the  footprint  looks  like  a  circular  sector. 


After  validating  the  6-DOF  model,  Clothier  then  works  on  a  footprint  approximation  technique. 
Clothier  uses  the  two  half  ellipses  to  model  the  basic  shape  of  a  footprint.  The  model  is  then  used 
to  approximate  the  extremities  of  the  true  footprint.  This  model  is  illustrated  in  Figure  1 
(reference  1). 


x 


Figure  1 :  Illustration  of  dual  half-ellipse  geometric  footprint  primitive.  Note  that  A  is  the 
initial  aircraft  location  at  the  time  of  failure,  AC  is  the  turn  radius  rt  for  the  chosen  bank  angle, 
CD  =  al,  CB  =  bl  =  b2,  EC  =  a2.  C  is  the  centre  point  for  the  frontal  and  back  ellipse. 
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In  Figure  1,  A  is  the  origin;  E,  B,  and  D  are  extremities.  Thus,  a\,  bi,  and  a2  can  be  used  to 
approximate  the  extremities  of  the  true  footprint.  C  is  the  center  of  both  ellipses,  and  C  is 
distance  rt  from  A. 


a1  =  maximum  distance  —  rt 
b1  —  rt  +  distance  of  glide 
a2  —  rt  +  distance  of  glide 

To  solve  for  ai,  bi,  and  an,  consider  the  front  ellipse  (yellow  portion)  bl  can  be  estimated  with  a 
descent  trajectory  comprising  of  a  0  =  90°  turn  and  glide.  This  trajectory  is  illustrated  in  Figure 
2  (reference  1). 


Figure  2:  Illustration  of  geometric  approximation  of  descent  trajectory  showing  (a)  3D  view, 
(b)  top-down  (xy  plane)  view,  and  (c)  side-on  view.  This  trajectory  comprises  a  turn 
(modeled  with  a  circular  arc  AB  of  radius  rt),  and  a  straight  line  glide  BD.  A  constant  FPA 
is  assumed  for  the  turn  and  similarly  for  the  glide.  AE  and  BE  are  tangential  to  the  arc  AB. 
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Using  Figure  2  and  assuming  a  level  turn,  you  can  approximate  the  turn  radius: 

V2 

r  = - 

g  tan/r 

Using  Figure  2(b),  it  is  possible  to  determine  the  distance  traveled  in  the  turn: 

dt  =  0rt 

Using  Figure  2(c)  of  the  illustration,  it  is  possible  to  find  AF  =  5 h,  given  AB  =  d,  assuming  yt  is  a 
small  angle: 


Sht  =  dt  tanyt  *  dtyt 

Where  yg  is  the  FPA  in  the  glide  or  the  inverse  of  the  maximum  lift  to  drag  ratio,  FPA  in  the  turn 
(Yt)  is: 


Yt  ktgYg 

Going  back  to  the  original  assumptions,  ktg  =1.5  for  a  bank  angle  of  45  deg.  Constant  ktg  is 
determined  empirically  and  enables  accurate  estimates  of  bi  and  a2  which  is  shown  below: 

H  —  Sht  H  -  Sht 
dg  ' 

tan  Yg  Yg 

bi=rt  +  dg 

a2  —  rt  +  dg 

al  —  djy[ax  U 

This  approximation  only  works  for  altitudes  where  a  90  deg  turn  is  possible.  Another  method  for 
finding  bi  is  needed  when  dg  <  0.  bi  is  an  approximation  of  the  maximum  across  track 
displacement.  Using  simple  geometry,  across  track  displacement  due  to  the  glide  and  due  to  the 
turn,  for  0  <  0  <  90°: 


8yg  =  dg  sin  0 

0 

8yt  =  dt  sin  — 

0 

dt  —  2r  sin  — 
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This  means: 


dg 


H  -  2rsin|tanyt 


tan 


Yg 


0 

8y  =  dq  sin  0  +  dt  sin  — 


8y  — 


tf-2rsin|tanyt 
tan  Yg 


+  2r  sin2 


0 

2 


H  2rtanyt  0  0 

5v  = - sin  0 - sin  —  sin  0  +  2r  sin2  — 

tan  yg  tan  yg  2  2 


H  2  T* 

Let  K\  = - ,  /G  = - and  K-,  =  2 r,  to  detennine  displacement  find  the  differential  of 

1  tan  Yg  2  tanyfl  3  F 

8y  with  respect  to0.  Then  to  detennine  if  the  displacement  is  a  maximum  find  the  second 
derivative  of  8y  with  respect  to  0.  Clothier  does  an  example  of  this  and  proves  that  the  solution 
is  a  maximum.  Therefore: 

bt  =  8y 


Given  the  above  methods  for  finding  al,  bl,  and  a2  (Figure  1),  it  is  possible  to  find  the  total 
footprint  area  A: 


n  n 

A  =  -aA  +-a2b1 

Clothier  than  validated  empirically  his  approximation  against  the  simulation  results,  He  found 
that  his  method  produced  an  accurate  approximation  of  the  true  footprint  extents  (see  reference 
1)  for  the  details  on  the  simulation  validation). 

Clothier  also  found  that,  at  low  altitudes,  the  footprint  area  was  more  of  a  pie  slice  shape  or  a 
circular  sector.  This  is  due  to  the  loss  of  altitude  during  the  turn  prohibiting  the  UA  from 
reaching  some  areas  behind  it  at  the  time  of  failure.  The  approximation  for  this  pie  shape  crash 
location  is  provided  in  (reference  1). 

SENSIS  METHOD 


Another  method  of  finding  the  Crash  Location  Footprint  for  a  fixed  wing  aircraft  is  the  Sensis 
Method  (reference  2).  Instead  of  doing  six  differential  equations  to  form  a  geometric  shape,  the 
Sensis  group  determines  an  outer  edge  and  an  inner  edge  to  form  the  polygon.  The  outer  edge  of 
the  footprint  is  limited  by  the  maximum  glide-down  range  of  the  aircraft;  making  sure  to  take 
into  account  the  loss  of  altitude  due  to  a  heading  change.  The  inner  most  edge  of  the  polygon  is 
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restricted  by  the  turn  radius  of  the  vehicle,  which  is  a  function  of  the  bank  angle  and  speed.  An 
example  of  this  footprint  can  be  seen  in  Figure  3  (reference  2). 


Figure  3:  Sensis  Method  Crash  Location  Footprint 

Although  the  Sensis  method  is  another  way  of  defining  a  crash  footprint,  there  are  reasons  that  it 
is  not  a  good  method  for  this  effort.  Primarily,  this  method  focuses  on  the  UAS  operator  finding 
a  termination  point,  which  implies  that  the  failure  is  not  catastrophic.  That  reasoning  excludes 
the  area  immediately  underneath  the  aircraft  and  behind  the  aircraft,  which  are  logical  locations 
where  an  UA  could  crash  if  there  was  a  catastrophic  failure.  However,  the  Sensis  method  does 
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confirm  some  of  the  methodology  used  in  Clothier’s  model  in  that  the  outer  edge  of  the  footprint 
is  limited  by  the  maximum  glide-down  range  of  the  aircraft. 

HELICOPTER  MODEL 


A  search  of  published  literature  found  no  current  helicopter  model  that  can  find  the  crash 
location  using  flight  information  and  different  types  of  crash  trajectories.  The  two  main 
assumptions  that  can  be  made  for  a  helicopter  crash  is  that  there  is  no  autorotation  and  drag  is 
neglible1.  With  these  assumptions  the  helicopter  crash  trajectory  becomes  projectile  motion.  This 
makes  lift  and  drag  negligible  and  also  means  that  majority  of  the  force  of  the  crash  will  be 
straight  down.  Since  the  height  and  the  speed  of  the  UAS  are  already  known  the  location  of  the 
crash  can  be  detennined  using  the  following  equation. 

R  —  distance  from  Helicopter  loss  of  control  to  Helicopter  Crash  site 

m 

g  =  9.81—  ( gravity ) 

V0  —  Helicopter  Speed  at  loss  of  control 
h  —  Helicopter  Height  at  loss  of  control 


R  =  V0 

HELICOPTER  MODEL  WITH  AUTOROTATION 


If  the  UAS  has  autorotation,  this  model  changes  drastically.  The  researchers  investigated  a 
method  that  can  be  used  to  apply  autorotation  to  a  UAS  helicopter  crash  model  (reference  3). 
This  method  states  that  gliding  descent  for  helicopters  rely  on  the  lift  generated  by  the  rotors. 
Autorotation  is  the  means  of  an  unpowered  rotor  to  maintain  rotor  RPM,  lift,  and  control.  The 
helicopter  behaves  such  that  the  relative  wind  comes  upward  through  the  rotor  disk.  The  aircraft 
is  giving  up  altitude  (or  potential  energy)  at  a  controlled  rate  for  kinetic  energy  to  drive  the  rotor 
as  shown  in  Figure  4  (reference  4). 


1  In  reality,  the  tumbling  helicopter  would  have  very  high  and  non-constant  drag,  but  this  assumption  provides  the 
most  conservative  estimate  on  how  far  the  helicopter  could  go. 
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Figure  4:  Momentum  Theory  of  Autorotation  through  a  Disc;  Illustration  by  John  Ball. 

A  simple  way  to  model  the  autorotation  for  the  large  variety  of  shapes,  sizes,  and  configurations 
of  rotorcraft  is  with  a  momentum  model  (Figure  5,  reference  4).  The  momentum  model  is 
presented  by  Leishman,  Johnson,  and  Prouty  (references  4,  5,  and  6).  A  momentum  theory 
solution  states  that  the  following  assumptions  are  made:  for  axial  flight  (motion  only  in  the 
vertical),  the  rotor  disk  can  be  modeled  as  an  infinite  number  of  blades  and  that  the  disk  is  strong 
enough  to  support  pressure  differences.  As  the  rotorcraft  descends,  the  relative  flow  field 
velocity  far  below  the  aircraft  is  the  descent  velocity.  At  the  rotor  disk,  the  velocity  of  the  air  is 
reduced  by  the  downward  velocity.  This  downward  velocity  is  the  hover  induced  velocity.  This 
slows  the  airstream  as  it  moves  past  the  disk.  Downstream  of  the  rotor  (which  in  this  case  in  the 
up  direction),  the  air  separates  and  results  in  a  relative  velocity  of  the  descent  velocity  minus 
twice  the  hover  induced  velocity. 


!k  !k 

Figure  5:  Solutions  to  Momentum  Theory  Autorotation  and  Young's  Approximation 
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The  physical  solutions  of  this  problem  are  complex  and  can  be  seen  in  the  image  above.  Young 
approximated  the  ratio  of  descend  velocity  to  induced  hover  velocity.  Johnson  built  upon  this 
and  included  an  additional  term,  figure  of  merit.  The  result  is  a  straightforward  equation 
(reference  5). 


TJT/2pA 
FM  =  — 

'Hover  Actual 

Where  FM  represents  the  figure  of  merit  and  k  is  the  rotor  efficiency.  The  figure  of  merit,  as 
defined  by  Johnson  is  the  ideal  power  required  to  hover  over  the  actual  power  required.  In  a 
steady  descent,  the  thrust  can  be  approximated  by  the  weight  of  the  aircraft  and  the  area  as  the 
area  of  the  blades.  If  the  power  required  for  hover  is  not  known,  more  advanced  helicopter 
design  should  have  figures  of  merit  between  0.75  and  0.85  so  a  conservative  assumption  of  0.7  is 
used.  Analysis  shows  that  the  figure  of  merit  is  significantly  more  influential  in  the  value  of  the 
ratio  of  vertical  descent  velocity  to  hover  induced  velocity.  Therefore,  the  model  defines  k  —  1, 
or  an  ideal  rotor  efficiency. 


*D  axial 

Vi 


FM  K  _  7K 
1+3 K  1+3 K 


The  induced  hover  velocity  can  be  calculated  using  equation: 


This  process  can  help  estimate  the  descent  rate  of  a  helicopter  in  axial  flight  conditions.  Where  p 
is  the  air  density  (assume  constant  sea  level),  T  is  the  thrust,  and  A  is  the  blade  area. 

Once  the  velocity  of  descent  is  known  the  distance  traveled  can  easily  be  calculated  by  replacing 
gravity  with  the  decent  velocity  or  with  the  rate  of  decent  velocity. 


TARGET  LEVEL  OF  SAFETY  MODEL  TOOL 


For  the  TLS  model,  some  assumptions  were  made  to  simplify  the  calculation  methods  for  PCL 
so  that  the  calculations  and  the  output  could  be  performed  by  a  tool  built  to  use  as  few  variables 
as  possible.  The  Clothier  method  was  chosen  as  the  model  for  fixed  wing  aircraft.  The  tool  uses 
glide  ratio,  aircraft  max  speed,  and  aircraft  weight  to  do  the  necessary  calculations.  It  was 
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decided  for  the  TLS  tool  to  implement  both  the  ellipse  and  pie  slice  models  and  have  the  tool 
switch  between  the  ellipse  and  pie  slice  patterns2.  For  helicopters,  the  assumption  was  made  that 
autorotation  was  not  an  option  based  on  the  capabilities  of  current  unmanned  helicopters. 
Therefore,  the  tool  will  use  the  projectile  model  for  all  helicopter  failures. 

The  crash  location  calculations  are  just  one  portion  of  the  TLS  tool,  the  other  portion  is  LCA. 
Although  the  LCA  is  not  important  to  the  PCL  area,  and  therefore  not  described  in  this  paper,  the 
crash  location  area  does  pull  the  operating  velocity,  and  the  max  glide  ratio  from  that  portion  of 
the  tool.  Figure  6  is  a  photo  of  the  crash  location  area  page  in  the  TLS  tool. 


Select  Aircraft  TvDe  HUafiai^BS 

■ 

General  Info: 

h  DM 

65000  ft 

Q|,H 

■3 

112  Knots 

If  helicopter 

Rotor  radius  DJ 

1 

12.5  ft 

If  Airplane 

Max(L/D)QH| 

►  I 

11.10374 

Distance  from  Helicopter  LOC  N/A  miles 

Area  of  Heliocter  Crash  N/A  miles2 

Area  of  Airplane  Crash  58244.476  miles2 

at=  136.48361  miles 

bt=  136.40875  miles 

a2=  135.91352  miles 


Figure  6:  TLS  Tool  Crash  Location  Page 


2  The  switching  altitude  depends  on  several  aircraft  parameters  including  airspeed  and  lift  to  drag  ratio.  It  typically 
occurs  around  50  to  100  ft  AGL  altitude  although  very  fast  aircraft  maintain  the  pie  slice  model  up  to  altitudes  of 
~1,000  ft.  The  exact  algorithm  for  this  switching  is  documented  in  the  3PRAT  report. 
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